###################################################################
###################################################################
## Required libraries
###################################################################
###################################################################

rm(list=ls())
source("Functions_Criterion.R")


library(msm)
library(MASS)
library(truncnorm)
library(VGAM)

###################################################################
###################################################################
## Set parameters for saving and length of rn
###################################################################
###################################################################

##Set to TRUE to do a long run to replicate what's in the paper
longrun<-FALSE

##Set to TRUE to save the result from the loop
saveoutput<-FALSE

###################################################################
###################################################################
## Analyses
###################################################################
###################################################################

##Generate a set of statistics for a single run
c1<-crit.function(9,longrun)

##Generate a set of statistics for all runs (note: published results were parallelized on a cluster)

tab.run<-NULL
for(i in 1:96) tab.run<-rbind(tab.run,crit.function(i,longrun))

C<-tab.run

if(saveoutput) save(C,file="Criterion")